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Abstract 



Explicit exact formulas are presented, up to fourth order in a strict chi- 
ral covariant derivative expansion, for the normal parity component of the 
Euclidean effective action of even-dimensional Dirac fermions. The bosonic 
background fields considered are scalar, pseudo-scalar, vector and axial vector. 
No assumptions are made on the internal symmetry group and, in particular, 
the scalar and pseudo-scalar fields need not be on the chiral circle. 
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I. INTRODUCTION 



The effective action plays an important role in quantum field theory, both from the 
phenomenological and from the formal points of view, since it compactly embodies the 
renormalized properties of the system; Green's functions, 5*-matrix elements and the expec- 
tation values of observables can be extracted from it. In addition it is suitable to study 
non-perturbative issues. Consequently the effective action functional has been the subject 
of a very large amount of work from different points of view, which include its proper defi- 
nition and renormalization, its symmetries and their anomalous breaking, its properties in 
curved space-times or finite temperature, and its calculation. An excellent and quite com- 
plete review for chiral fermions in flat space-time and zero temperature, the case of interest 
to us, can be found in |1J]. For anomalies in curved space-times see e.g. 0. Work on chiral 
fermions at finite temperature can be found in 

In the present work and its companion paper |J we concentrate on the computational 
issues. Concretely we deal with the calculation of the effective action of Dirac fermions in the 
presence of bosonic external fields, thus the functional adds one-loop Feynman diagrams with 
fermions running on the loop and bosonic external legs. We restrict ourselves to the case of 
even dimensions (since there are some technical differences with the odd-dimensional case), 
zero temperature and flat space-time. The class of external bosonic fields to be included is 
that of scalar, pseudo-scalar, vector and axial vector fields. Coupling to higher tensor fields 
is not included, nevertheless, the class considered here is quite large since no assumption 
will be made on the internal symmetry group, that is, the external fields are arbitrary 
matrices free from algebraic assumptions regarding their dimension, commutativity, chiral 
circle constraint or any other constraints among them. 

The computation of the effective action functional in closed form is not possible in general, 
and thus several asymptotic expansions have been devised. In the heat kernel expansion the 
terms are classified by its scale dimension, that is, each term has a well defined number 
of external fields and derivative operators. It is the computationally simplest expansion 
and so it has been carried out to considerably large orders even for curved space-times 
Other expansions can be regarded as resummations of this one. In the perturbative 
expansion, the contributions are classified by the number of external legs, that is, the number 
of external fields in the term, and all orders in their momenta are added. This is a weak field 
approximation which however captures the non-locality of the exact functional. On the other 
hand, in the covariant derivative expansion, to be considered in this work, the contributions 
are classified by the number of chiral covariant derivatives, or equivalently, by the number of 
Lorentz indices. Both definitions are equivalent in the absence of external tensor fields. This 
counting is appropriate for external fields with a smooth space-time dependence and weak 
gauge fields. The scalar and pseudo-scalar fields need not be weak, and this allows to study 
non-perturbative issues such as spontaneous symmetry breaking. Two important properties 
of this expansion is that chiral symmetry is preserved (modulo anomalies) separately for 
each term, and that the terms are local. 

The effective action is a ultraviolet divergent quantity which needs to be renormalized. 
The non-perturbative definition of this functional in the general chiral case is not so straight- 
forward as for vector-like theories (in which there is no pseudo-scalar nor axial fields) due 
to the presence of essential chiral anomalies which affect the imaginary part (the phase of 
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the fermionic determinant) |7j]. Nevertheless such a definition exists and so the effective 
action is perfectly well-defined also in the chiral case, displaying only the standard renormal- 
ization ambiguities in the form of polynomial terms. The subtleties in the non-perturbative 
definition of the effective action are much alleviated within the asymptotic expansions noted 
above. In particular, within the covariant derivative expansion to be worked out in this 
paper for the real part of the effective action and in || for the imaginary part, there is only 
a finite number of ultraviolet divergent terms. They are afflicted by polynomial ambiguities 
(including anomalies and multivaluation) but are otherwise unique. All higher order terms 
are ultraviolet finite and so unambiguous. 

It should be emphasized that the covariant derivative expansion of the effective action 
functional, although asymptotic, is perfectly well defined, that is, it does not depend on how 
it is written or computed. (This is not true for other expansions such as the commutator 
expansion, as we will discuss below.) This is because it corresponds to classify the terms 
by their scaling under (covariant) dilatations of the external fields. Each of these terms is a 
well defined and universal functional, in the sense that they hold for any possible internal 
symmetry group, and the same is true for the perturbative or heat kernel expansions. For 
these two latter expansions it is relatively easy to write an explicit form for each term without 
putting restrictions on the internal symmetry group. The reason for this is that all non 
commutative quantities, namely the matrix- valued external fields, are treated perturbatively 
and each term contains only a finite number of them. On the other hand, in a strict covariant 
derivative expansion the gauge fields are treated perturbatively but the scalar and pseudo- 
scalar fields appear to all orders in every single term, and these two fields do not commute 
with each other. Nevertheless, in this work and in || we show that the terms of the derivative 
expansions are also fully amenable to explicit computation in closed form without assuming 
particular properties of the internal symmetry group. In addition, no chiral rotation (or 
diagonalization) is required to express these universal functionals. Our result takes an 
analytical form in terms of the external fields. 

In this work we deal with the real part of the effective action and compute it until fourth 
order in the covariant derivative expansion for arbitrary even space-time dimensions. The 
real part is the simplest one since it is free from essential chiral anomalies and multivaluation. 
It contains only a scale anomaly. In || the study of the imaginary part is dealt with at 
leading order in the derivative expansion for two- and four-dimensional space-times. Besides 
the results themselves, in this paper and in || we introduce notational conventions which 
are very well suited to the chiral problem. In particular, we find that the formulas for the 
real part in the full chiral case are identical to those of the vector-like case. 

In Section [TJ we introduce our notational conventions, some of which are not standard, 
and also introduce the effective action. In Section |TTl| we show how, for the real part of the 
effective action, an appropriate notation allows to reduce the full chiral case to the vector-like 
case, and also to carry out explicitly all integrations over the momentum in the fermionic 
loop. The Section ends with explicit formulas for the real part of the effective action up to 
four covariant derivatives, based in the convenient method introduce by Chan for bosons 
H . Finally in Section [IV| we illustrate the meaning of the formulas by analyzing the case of 
second order and two space-time dimensions and a particular case is worked out explicitly. 
Next we show how our notation allows to obtain commutator expansions quite efficiently, 
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and the Section is ended by giving the analogous explicit formulas for bosons, once again 
without restriction on the internal symmetry group. 

II. GENERAL CONSIDERATIONS 

A. The Dirac operator 

The Euclidean effective action of fermions in a (i-dimensional flat space (d even) is / d d xipT)il), 
where D is the Dirac operator. The class of operators to be considered is, in terms of the 
left-right (LR) fields, 

D =I) R P R + I) L P L + m LR P R + m RL P L , (1) 

where P Rt L = |(1 ± 75) are the projectors on the subspaces 75 = ±1. Our conventions are 

7 M = 7j > {Ta» lu} = 28^ , 75 = 7] = % l = r/ d 7o ■ • ■ 7d-i , tr Dirac (l) = 2 d/2 . (2) 

With r/d = ±i d / 2 (a concrete choice will not be needed). D^' L = <9 M + vj^' L are the chiral co- 
variant derivatives. The external bosonic fields v^ ,L (x) and rriL R (x), mm(x) are matrices in 
internal space (referred to as flavor), the identity in Dirac space and multiplicative operators 
in x space. Unitarity of the theory imposes restrictions on the hermiticity properties of these 
fields, namely, v r> l must be antihermitian and m) LR = mm- In practice these restrictions 
will play almost no role in the calculation^ and and m R L will be treated as independent 
variables. To avoid infrared divergences we will assume that vtilr and mm are nowhere 
singular matrices. 

In terms of the vector-axial (VA) variables 

D=^ v + 4j 5 + S + j 5 P, (3) 

where DY = + is the vector covariant derivative and 

v R , L = V±A, m LR = S + P, m RL = S-P. (4) 



B. Chiral transformations 

Chiral transformations act as follows 

m LR -> Vtl l m LR Vt R , m RL -> Q^mmQi , (5) 



1 However we will exploit the fact that the spectrum of m R (x) (defined below in eq. (g)) is non 
negative to choose the branch cuts of the logarithm or squared root functions along the negative 
real axis. 
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where Q R} l{x) are independent, nowhere singular and otherwise arbitrary matrices in flavor 
space. Vector gauge transformations correspond to chiral transformations in the diagonal 
subgroup fl R = Ql, i.e. D — > fi _1 DQ. 

It will be convenient to introduce the two combinations 

m R = m RL m LR , rn\ = m LR m RL , (6) 

which transform solely under Q R or respectively. Note that the two matrices m RL (x) 
are related by a similarity transformations and thus they have the same spectrum. 

The different pieces in the Dirac operator transform in a well defined manner under 
chiral transformations, namely, the quantities ttilr, m R L, D R and Dl fall in the chiral 
representations LR, RL, RR and LL respectively, cf. eq. (||). New objects with well 
defined chirality are obtained by multiplication in the natural way, i.e., if X ab falls in the 
representation ab and Y bc in be, for a,b,c = R, L, the product Z ac = X ab Y bc (no sum over b 
is implied) falls in the representation ac. If in addition X ab is a multiplicative operator, its 
chiral covariant derivative 

(D li X) ab = D*X ab -X ab B b fl , a,b = R,L (7) 

is also multiplicative. In particular, 

Df.rriLR := (D^m) LR = <9 M m L/? + v^m LR - m LR v^ , 

D^m RL := (D^rl = d^rriRL + v*m RL - m RL v^ . (8) 

In addition, 

F* = iD«,D*], F^ v — [D^, D^] , (9) 

are also chiral covariant and multiplicative. All multiplicative chiral covariant local objects 
come as combinations of m, F and their chiral covariant derivatives. 

C. The effective action 

The fermionic effective action W = — log/ VipVip exp(— J d d xip~Dip) is given by 

W[v,m) = -Tr log(D). (10) 

As is well-known this expression is formal due to the presence of ultraviolet divergences. 
Mathematically proper definitions of Tr log(D) exist in the literature. Here it is only nec- 
essary to emphasize that different definitions of W may at most differ by terms which are 
local polynomial of dimension d, that is, polynomials in the external fields v and m and 
their derivatives. This is a standard result of perturbative quantum field theory that can 
be established by isolating the ultraviolet divergent one-loop Feynman graphs, which may 
contain at most d insertions, and expanding on the external momenta to extract the diver- 
gent part. In practice this means that any method consistent with the formal expression 
can be used to make W finite, since it will give the same ultraviolet finite contributions as 
any other method. The effective action of a concrete physical system described by the Dirac 
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operator D will be given by any of the renormalized versions of W plus an appropriate local 
polynomial counterterm. 

As usual, it will be convenient to introduce the pseudo-parity transformation, R <-> L 
(that is, vl <-> v R and m LR <-> m RL or equivalently A M — > —A^, P — > —P) and split the 
effective action into its even and odd components under this transformation, W^, i.e., 

W[v,m] = W + [v,m] + W~[v,m] . (11) 

(Also known as normal and abnormal parity components.) Due to parity invariance, which 
involves and additional (xq, x) — > (xq, —x) in the fields, the pseudo-parity odd component of 
the effective action W^fi^m.] is that containing the Levi-Civita pseudo-tensor. In addition, 
W^l/u, m] coincide with the real and imaginary parts, respectively, of the (Euclidean) effective 
action when the standard hermiticity for the fields is assumed. 

By computing the effective action within a derivative expansion, we mean expressing 
it as a sum of terms with a well-defined number of covariant derivative operators but any 
number of scalars. This is equivalent to classify the terms by the number of Lorentz indices 
they contain. We emphasize this since sometimes the expression "derivative expansion" is 
used in the literature to denote large mass expansions or expansions in the number of all 
kind of commutators, that is, those implied by the covariant derivative plus those of the 
form [m 2 , ], etc. Further, we mean writing each of the terms using building blocks which 
are multiplicative operators in x space (i.e., functions of x and not differential operators) 
and trivial in Dirac space (since the Dirac trace can be explicitly computed for any given 
order in a derivative expansion). The trace in flavor space will not be worked out since we 
are not imposing any algebraic constraint on the flavor structure of the fields. 

Note that, because d is even and all Lorentz invariants are formed with 5^ and eox...d-i, 
there are no odd-order terms in the derivative expansion. 

In this work we will consider the pseudo-parity even component of the effective action 
which is simpler due to its lack of chiral anomalies. The pseudo-parity odd component is 
worked out in ||. 



D. Local basis in flavor space 

Occasionally it will be useful to diagonalize the flavor matrices rriL R (x) and m R L(x). 
This can be done by solving the eigenvalue problem (at each point x) 

which yields 

m LR \j, R) = rrijlj, L) , m RL \j, L) = m^j, R) , 

(j, L\m LR = rrijij, R\ , (j, R\m RL = m^j, L\ . (13) 

{j, R\ is the dual basis of \j, R), (j, R\k, R) = 5jk (no orthonormality of the basis is implied). 
The numbers m 2 - are the common eigenvalues of m 2 R and m| and \j, R) and \j, L) are their 
eigenvectors. Because rn) LR = m R L, the eigenvalues rrij can be taken to be positive. 
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E. Some notational conventions 

In this section we will introduce some notational conventions which are essential for 
carrying out the subsequent calculations. They are also used in 0. 

Because W + is by definition invariant under the exchange of the labels R and L, each 
term T in the expansion of W + will have a pseudo-parity conjugate term, T*, obtained from 
T by exchanging everywhere the labels R and L. (Note that due to the cyclic property, it 
may actually happen that T and T* coincide.) Thus we will use the following convention: 

Convention 1. In W + , the terms T and T* will be identified, so that under this Convention 
T actually stands for |(T + T*). 

It is important to note that word "term" is used here in a very specific way, namely, it refers 
only to contributions of the form T = tr(X), that added produce the effective action. (By 
extension, term may denote also the quantity X itself. In what follows two such quantities 
X and Y differing only by a cyclic permutation will be identified, since they are equivalent 
inside the trace.) Of course, the identification implied in Convention 1 does not apply to 
smaller pieces or factors inside each term. 

Due to chiral covariance each term is a product of factors with well-defined chirality 
(namely, LR, RL, RR or LL) correctly combined to preserve chirality (i.e., • ■ -X^Y^- ■ ■). 
In addition, due to the cyclic property of the trace, if a term starts with label a — R, L it 
has to end also with same label a, e.g. tr(X rr Y R lZlr). Based on these observations we 
will make the following convention: 

Convention 2. In expressions where the chiral labels are combined preserving chirality, 
these labels are redundant and will be suppressed, thus a term such as ti{X RR Y RL Z LR ) 
will be written as tr (XVZ) . Note that tr(XVZ) could be expanded either as ti{X RR Y RL Z LR ) 
or tr(X llYlrZ R l), but both expressions are equivalent under the Convention 1. We 
will choose the first label as R. 

For instance 

tr(F> D^m D v m) = tr(F^r^m KL A,m Li? ) 

= -tT{F^b^m RL b v m LR ) + ^tr(F^ry™ LR A,m RL ) . (14) 

This notation allows to use the objects m, etc, as ordinary operators (i.e., elements of 
an associative algebra). For instance, the property mf(m) = f(m)m is verified as is readily 
checked. With this notation eqs. (||) and (|9]) become 

D M m = D^m - mD^ = [D M , m] , F ^ = [D„, D v ] , (15) 

and eq. ( |l3l) becomes 

m\j) = m j \j), (J\ m = m j (J\. (16) 
Another essential property, the cyclic property of the trace, holds too. For instance, 

tr(Xm) = tr(X RL m LR ) = tr(m LR X RL ) = tr(m RL X LR ) = tr(mX) . (17) 
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Let us note that these conventions have to be slightly modified to include the pseudo-parity 
odd component W~[m, u]. In particular, the cyclic property of the trace is modified and 
this fact is at the origin of the chiral anomaly in this formalism 0. 

Below, it will be necessary to carry out parametric integrations where the parameters 
appear in different places tied to operators which do not commute. To do these integrations 
another convention is needed regarding the order of operators. Consider for instance the 
integral 



d d p 

tr 



(2vr 



m m 



p 2 + m 2 11 p 2 + m 2 M 



(18) 



It can be computed by first taking the trace using the basis of eigenvectors of m introduced 
in subsection [11 D| . This yields an expression involving matrix elements of in these basis, 
namely, 



d d p x , rrij m k 
(2tt)<* 2 1 <p 2 + m 2 ^ V,x)jk P 2 +ml [V ' J ^ )kj, 



(19) 

, j— / ( i — i j — — — I I ) ~ 

3,k 



where (v^jk = (j, L\v^\k, L) and (v^kj = (k, R\vjf\j, R). (For definiteness we have assumed 
that the operator inside the trace in eq. ( |TBD has been expanded as a RR term.) The point 
of taking matrix elements is that, since now everything is in terms of commuting numbers, 
the momentum integral can be done explicitly. Instead of that, we will use an equivalent but 
preferable procedure, namely, we will label the matrices m according to their position relative 
to v^. Since there are two i/ M there are three possible relative positions which are labeled 
by 1, 2 and 3. With this prescription, the same integral can be represented unambiguously 
by the formula 



tr 



(2tt 



mi m 2 



2 

2 ^2 , m 2 V » 



p A + mf p z + m 



(20) 



The usefulness of this notation is that, since the ordering is given by the labels, m\ and vrt2 are 
effectively c-numbers and the momentum integration can be carried out straightforwardly. 

Convention 3. In an expression f(Ai, B2, . . .)XY ■ ■ ■ the ordering labels 1,2,... will denote 
the actual position of the operators A,B,... relative the fixed elements X,Y, . . . so 
that A is to be placed before X, B between X and Y, etc. That is, for a separable 
function f(a, b, . . .) = a(a)/3(b) ■ ■ •, the expression stands for a(A)X(3(B)Y ■ ■ ■ 

Note that the fixed elements X,Y, . . . appear only as simple factors (i.e. perturbatively) 
whereas the labeled operators A,B,... may appear with a non perturbative functional de- 
pendence. Also, note that several different labeled operators can have the same label pro- 
vided that they commute. f\ Finally, in an expression with n fixed elements inside the trace, 



2 The separation into fixed elements and labeled operators is a matter of convenience. It is also 
possible to label all operators to indicate their relative position and this is often useful in order to 
carry out algebraic manipulations with a computer. 
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the labels 1 and n + 1 are equivalent, due to the cyclic property of the trace. So, in the 
previous example, m 3 = mi. We remark that Convention 3 is independent of Conventions 1 
and 2. 

The usefulness of this notation can be further exposed through the following observation 
regarding commutators. Let X be a single fixed element and let the symbol Da denote the 
operation [A, ], thus it immediately follows from Convention 3 that Da = A\ — A2, that is 

f(D A )X = f(Ai-A 2 )X. (21) 

For instance, 

e A Xe~ A = e M e~ M X = e M ~ M X = e° A X , (22) 

which is a well-known identity. (Alternatively, this example can be regarded as a proof 
of eq. (pip.) This observation will prove useful to carry out commutator expansions (cf. 
subsection II V B|) . 



As another illustration of Convention 3, consider the following identity, where f(A) 
depends on the operator A and 5 A represent some first order variation of it, 

6f(A) = f{Al J- f } A2) 5A. (23) 

Ai — A 2 

Here 6 A is the fixed element referred to in the Convention and Ai and A 2 refer to A before 
and after 5A, respectively. This identity can be proven as follows.0 As is well-known 

Se sA = f dte tA 5Ae {s ~ t)A , (24) 
J 

where s is a c-number parameter. Using Convention 3, this can be rewritten as 

5e sA = / dte tM+ ^ s ~ t)M 5A = — — 5 A , (25) 

Jo A x -A 2 

where in the last step we have used that A x and A 2 are commuting quantities. This identity is 
then generalized for arbitrary f(A) by Fourier transform. From another point of view, noting 
that Ai—A 2 is equivalent to [A, } for any A, eq. ( PB"D is equivalent to [A, Sf(A)) = [f(A), 5A], 
which is a trivial consequence of S[A, f(A)] = 0. 

Inside a trace the two ordering labels 1 and 2 become identical due to the cyclic property, 
therefore eq. (|23| ) yields 

Tr(5f(A))=Tr(f'(A)5A) . (26) 

(/' denotes the derivative of /.) The corresponding formula for two successive first order 
variations, is 

Tr (5'5f(A)) = Tr \f{A)5'8A + ^^"{^^ 'm) , (27) 
V A\ — A 2 ) 

which in particular implies Tr ([5', 5}f(A)) = Tr (f'(A)[S', 5]A). 



3 Alternatively it can be proven by starting from 6(A^ 1 ) = — A~ l 5AA~ l = —A± 1 A 2 1 6A and 
applying it to f(A) = (2-Ki)~ l J r dzf(z)/(z — A) where T is positively oriented and encloses the 
spectrum of A, or also it can be established by considering functions of the form f(x) = x n . 
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III. CALCULATION OF W + TO FOURTH ORDER 



A. Reduction to a vector-like theory 

A first benefit of Conventions 1 and 2 is that the functional VF + [i>,m] will be formally 
identical to the effective action of a vector-like theory. That is, if Wv[V, 5] denotes the 
effective action functional when = P = 0, 

W + [v,m\ = W y [v,m], (28) 

provided that jy + [t>,m] and Wv[V, 5] have been renormalized preserving chiral and vector 
gauge invariances respectively. This is because when the functional W + , expressed in our 
notation in terms of v and m, is particularized to the vector-like case (i.e., vr = Vl, m LR = 
mm) no simplification occurs; it remains unchanged in our notation. (W~ vanishes for 
vector-like configurations.) As a consequence, W + for the general case can be obtained by 
computing only the effective action for vector-like configurations. 

It is an essential point of this discussion that we are considering field configurations 
which are free of constraints in flavor space: the calculation of W /+ [f,m] through W / v[y,S'] 
requires and S to be generic in flavor space for the formal replacements V — ► v and 
S — > m in the formulas to be well-defined. For instance, if the calculation of Wv[V,S'] 
were carried out in the particular case of commuting and S, [V^, S] would be identified 
with zero and this would result in an ambiguity in iy + [t>,m] by terms m], however 
[ v n, m ]nL = [V/j,,mRi] + {A^^rl}, which does not vanish even for Abelian flavor groups, 
due to the axial term. Therefore the proper procedure is first to compute H^vt^j m general, 
then make the replacements V — > v and S — > m, and afterwards particularize the resulting 
formulas to the case at hand. This replacement is a kind of analytical continuation from the 
vector-like to the chiral case. 



B. Chan's Method 

Here we present the most efficient method to compute W + . Further heuristic consider- 



ations are made in Section |V[ We have already reduced the problem to a vector-like one, 
with effective operator D +m. Next, the theory is reduced to a bosonic one using D^D 
as operator. This procedure is standard. 

D=[?+m, D t = 7 5 D 75 = -/;+m. (29) 

The second equality implies that D and are related by a similarity transformation and 
so they can be assigned the same determinant, thus 



W + \v,m] 



iTrlog(DtD). (30) 



Two remarks are in order. First, the notation assumes the fields have the standard 
hermiticity, but only the fact that is similar to D is essential. Second, within a given 
particular renormalization prescription (e.g. C-fimction) eq. (|30D will be correct modulo local 
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polynomial counterterms. This is entirely sufficient since, as noted above, the action only 
determines the effective action modulo a local polynomial of dimension d. 
The bosonic operator D^D can be worked out to give 

D + D = -D\ - - [0, m] + m 2 , 

= P 2 + U, (31) 

where 7^ = 5^ + 07^, and 

P M = D, , P 2 = -Pi , U = m 2 - ^b„m - -a^F^ . (32) 

(Note that P 2 is unrelated to the pseudo-scalar field P = (tiilr — mm) 1 2.) 

A manifestly gauge invariant form for the derivative expansion of the effective action of 
a general bosonic theory with Klein-Gordon operator of the form P 2 + U has been given 
in m to fourth order and in || to sixth order. Those formulas assume only that P M is a 
covariant derivative, i.e., of the form d^ + W^x), and U{x) a Lorentz-scalar field. The fields 
Wfj,(x) and U(x) are matrices in some internal space. To fourth order the result isfj 



Tr log(P 2 + U) = J ^#tr[ - log(iV) + P -N 2 

'-2Nl + (Ay\y 2 + 2{NN m f + iNF^NNpNv + (F^N 2 ) 2 

(33) 



d(d 
+ ••• 



where the following quantities have been defined 

N =-T^ji N^iP^N], N m = [P ll ,[P„N]], F, U = [P„P U \. (34) 

In addition, N* stands for (A" 2 ) 2 , and (N^Ny) 2 = N^N U N^N U , etc. The symbol tr includes 
all internal degrees of freedom. The dots refer to terms with six or more Pp. All terms are 
multiplicative operators (i.e. just functions of x) and manifestly gauge invariant. Because 
all terms are multiplicative, the trace cyclic property applies without restrictions. In what 
follows, terms related by a cyclic permutation will be identified when they appear inside the 
trace. 

Eq. ( p3[ ) can straightforwardly be applied to compute W + with the identifications 

N^DpN, N IHl = blN, Fp V = Fp V , (35) 

and A" and U given in eq. (|3"2"D. It is just necessary to note that in the bosonic case U is 
assumed to be of zeroth order, however, in the fermionic case U contains terms of first and 



4 A misprint in the formula presented in || was corrected in || . 
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second order in D M , so N has to be reexpanded in order to have a well-defined derivative 
expansion of W + : 

N- 1 = A + [/« + , 

A=p 2 + m\ U^^-j^D^m, = ~a /lv F ftu . (36) 

Before proceeding, two comments can be made. First, DD^ can be used instead of D^D. 
It corresponds to the replacement m — > — m in the previous expressions, so W + is an even 
functional of m. This is readily verified in the final expression. Second, D 2 would also 
produce another acceptable bosonic operator P' 2 + U', with P^ = D M + 7 M m. However, 
this form is not suitable for a derivative expansion because 7 M m is of zeroth order. This 
implies that low orders in the derivative expansion of the fermionic problem would pick up 
contributions from all orders in the bosonic derivative expansion. 

Combining the previous equations, and after taking the Dirac trace, the following result 
is obtained:^ 




2 d l 2 is the dimension of Dirac space and tr stands for trace on flavor space only. Two further 
notational conventions have been used. First, if I is an ordered set of Lorentz indices X^i 
denotes^] D^Xi, thus, for instance, 



5 We will use the notations W^~[t>,m] and W7j[f,ffi] where n specifies the order in the derivative 
expansion and d the space-time dimension. 

6 Note that in || this is denoted by Xj^. 
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= D^m , = D M (D v m) . (40) 

Second, the bosonic determinant and all manipulations used are invariant under the trans- 
formation defined by 

m.W = /n, D^D\ = -D^ (XYy = Y^X\ (41) 

(this is just the Hermitian conjugation when the fields have the standard hermiticity) . So, 
for instance, 

m i = > m U = m w > F U =- f ^- ( 42 ) 

Therefore, all terms in the previous formulas either are self-conjugated or come in conjugated 
pairs. In writing the formula for W/jv,???,] each term has been identified with its conjugate 
PJTHfl. For instance A~ 2 (m 2 ) M A _1 m^A _2 m Ml/ , actually stands for 

1 1 / 1 1 11 1 / 2^ 1 

(Recall that terms related by a cyclic permutation are identified due to the trace.) 
In deriving the formulas for W + the identity 



1 \ 1, 2 , 1 

A) =-A (m >"A 



(44) 



has been employed. In general it will not be practical to further expand (m 2 )^ as mm^ + m^m, 
since m 2 is a chirality preserving quantity and so it is better behaved than m. From the 
method used it follows that W + can be written in such a way that single factors m appear 
only through derivatives, i.e. m M , m^, etc, and underivated m are always squared. 



C. Momentum integrals 

Using the Convention 3, all momentum integrals in eqs. (|3~7| , |3~8| , |3~9l) can explicitly be 
computed. The ultraviolet divergences can be dealt with by using dimensional regularization. 
The basic integral is of the form 

r , ^ (47r) d "/ 2 r(rf72) f d d p ( 2V » 1 

W*, e) = iA7T)d/2T{d/2) Jj^-Ap) It ^ (45) 

where 

d = d-2e, k = ^ + r - 1 , A, t = p 2 + z 2 . (46) 

Here, zi, . . . , z n are independent non vanishing variables (to be replaced later by m). The 
integral can be done straightforwardly and gives 

/ - nk-e+l)T(-k + e) » k 

(4vr) d / 2 r(rf/2) £f J n,i ' 1 j 
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with 



R n ,i(zi, ...,z n ) = J[ 2 1 2 . (48) 



jjii Z j Z i 



Finally, expanding in powers of e yields 

1 /„ / I 



L n,kA 



(47r) d / 2 r(d/2) 



e 



H n , k +[-- log(^) Q n , fc + 0(e) . (49) 



where the following functions have been introduced 

n 

H n!k ( Zl , ji) = (-l) fc+1 log(^//i 2 )i4, * , 

i=l 

n 

Qn,k{ Z h ...,Z n ) = {-If J2( Z f) kR n,i ■ (50) 



8=1 



Several comments are in order here. 

(i) As indicated by the notation the functions H n k and Q Ut k do not directly depend on 
the space-time dimension. 

(ii) Following the standard practice, the parameter /i 2 has been introduced for conve- 
nience in the dimensional counting. The integral I n ,k,d does not depend on /i 2 and so its 
dependence cancels in the r.h.s. of eq. (f4"9"p. Under a scale transformation, R n i and Q n ,k 
transform homogeneously, but 



\z1 , H nM » X k+1 - n {H n , k - log(A)Q n , fc ) . (51) 



Therefore, Q n ,k represents the contribution of the integral to the scale anomaly. 

(hi) Whenever the integral I n ,k,d is infrared and ultraviolet convergent, the function Q n ,k 
must vanish, 

Qn,k = when 1 < k + 1 < n . (52) 

This is because as e goes to zero I n ,k,d has to remain finite. In addition, log(/i 2 ) no longer 
appears and there is no scale anomaly. When the integral is infrared finite (k > 0, which 
is the case of interest) but not necessarily ultraviolet finite, the functions Q n ,k are just 
homogeneous polynomials in the squared masses, zf, since taking a sufficient number of 
derivatives of I n ,k,d with respect to the any of the zj, turns the integral into a ultraviolet 
convergent one. 

(iv) The ultraviolet finite prefactor [4^/2^/2) ^ as been introduced in the definition of 
In,k,d for convenience, in order to remove renormalization dependent constant terms (the 
usual log(47r) — 7 of the MS scheme) in eq. (|49|). Those constant terms would come with 
Qn t k and would amount to further polynomial contributions to the effective action. 

(v) From the definition of I n ,k,di it follows that H n k and Q n ,k are completely symmetric 
functions of the m 2 , and furthermore, they remain finite as two or more of the squared 
masses become equal. {Q n ,k is a polynomial so the finiteness is obvious in this case.) 

(vi) The identity 
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11 1 / 1 1 

a; a; - _ ' zi-z) \k~ a~ 



(53) 



gives rise to the following recurrence relation 



T I \ _ In-l,k,d( z 2, Z3, . . . , Z n ] e) — I n -l,k,d(zi, %3, ■ ■ ■ , Z n \ e ) / KA \ 

J-n,k,d\Zi, Z2, ■ ■ ■ ; 6J — 2 2 ' l" 34 / 

z x — 2 2 

As a consequence, analogous recurrence relations apply to H n ^ and Q n ,k as well, and they 
can be used to compute these functions starting from 

H Xik = (-l) fc+1 (^) fc log(z 2 //i 2 ), Q x , k = (-l)W- (55) 

More generally, it will be necessary to consider also the case of A; raised to different 
powers in the momentum integral, eq. (flop, namely 



4, M K*x, ■■■,z n ,e)- (4 w /2r(d/2) _/ ^ J 11 A? ■ (5b) 



(47r) d / 2 r(d/2) J (2 



The Ti are assumed to be positive integers. Obviously, this integral is given by the previous 
formulas by taking the first r\ arguments to be all of them equal Z\, then the next r 2 
arguments to be z%, and so on. That is, 

In,k,d "(^1) ■ ■ ■ i z n ;e) = J nH \-r n> k,d{zi, ■ ■ ■ , Z\, . . . , z n , . . . , z n ; e) , (57) 

where in the r.h.s. each Zi appears times. The corresponding analogous definitions apply 

to H r n X' rn and g- r r "" 

Alternatively, the same result is obtained taking derivatives with respect to the zf. For 
instance, 



Hn,k" ,rn ( z l> ■ ■ ■ ) Zn, — 11 



(r i -l)!^?) ri - 1 



E n)k {z x ,...,z n \[i). (58) 



Analogous formulas hold for Z^ 1 /!"/™ and "' r ™. 

It is interesting to note that has the form of a rational function times a logarithm of 
zf I fj, 2 (see the remark (iv) above). This needs not be true for if^'"' r "; taking the coincidence 
limit Zi = zj or derivating with respect to zf yields further terms which are purely rational 
functions of the masses, without logarithms. 



D. Normal parity effective action through fourth order 

The momentum integrals in eqs. (|37| , |38|j39D can straightforwardly be worked out using 
the previous formulas for I^ k "/ n - I n the spirit of dimensional regularization, the variable 
d in eqs. (^,^,^) (excepting that in d d x) should be replaced by d = d — 2e and a Lau- 
rent expansion in e is to be performed. After that there will be two types of terms. First, 
contributions coming from (1/e — log(fi 2 ))Q r ^' k "' r ™ in the momentum integral. These pick 
up ultraviolet divergent terms plus finite terms from the e in d. And second, finite terms 
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from H^ k "' ,rn . The first kind of contributions are just polynomials and so renormalization 
dependent. Because all our formulas are already manifestly gauge invariant before momen- 
tum integration, those polynomials are gauge invariant and removable by counter terms. 
(Within other schemes, the divergent polynomial terms would still be gauge invariant but 
not necessarily the finite parts.) Therefore modulo polynomials, the effective action can be 
obtained by using d instead of d and keeping only the finite part H^ 1 k " ,rn in the momentum 
integrals. 

The calculation of Wo"[v, m] in eq. ( j3?|) can be done by first taking the derivative with 
respect to m 2 . This removes the logarithm and allows the application of the momentum 
integrals derived in the previous subsection. Wo~[i>, m] is then obtained by integrating back 
with respect to m 2 . This introduces an ambiguity which is just a polynomial. 



W+[v,m] 



_^d/ 2 2 d / 2 



2(47r) d / 2 r(d/2 + 1) 



d d xti 



m d \og{m 2 /n 2 



(59) 



Using our Convention 3 and the integrals of the previous section, the expressions for 
iy 2 + [' y ; m ] an d W 4 + [t>,m] are a direct transcription of eqs. ( |38l , |39"l) . In the following formulas 



the arguments z±, . . . , z n of H^ k "' n have to be substituted by mi, ... , m n . 



W£ d [v,m] 



1 



0<2/2 



2 (4vr) rf / 2 r(rf/2) 



d d x tr 



r 2.2 



2 H 2,d/2~i m l + \ H %4/2K m >V 



,2\2 



(60) 



W+ d [v,m] 



1 2 d / 2 
\(A<n) d / 2 T(d/2) 



d d x tr 



i r_rl,l,l r 
1 rr2,2 2 



+ \H#j£ 1 (m lt m v )' 
■ 4- - H 1 ' 1 F 2 

[iv ' A J1 2,d/2-l r nv 

j H 3 2 dM m2 )» m » m ^ 



tj-2,2,2,2 / 2\2/„,2\2 2 rj2,2,2,2 // 

d(dTT) H ^ {m } > ] » ~ d(d + 2) A > d,2+l{{ 

16 rr 3.1.3.1 /-.__2\2/.__2\2 j rr3,3 / 

Jy d(d + 2) 2 ^ 2+l1 



16 rr3,l,3,l / 2^ 

d(d + 2) 4 < d / 2+l1 ' 

i 16 rT3,l,3 / 2N2 

+ d(d+2) H ^ 2+l[m ^ 



8 2 2 2 
"rf(rf + 2) 3 ' rf / 2+1 



m 2 ),(m 2 



2 22 2 
d(d + 2) H2 ' d / 2+1 7 



(61) 
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IV. DISCUSSION 



A. Expression in terms of matrix elements 



In order to analyze the previous formulas, let us consider in more detail the second order 
term in two space-time dimensions, W / 2 + 2[ t '; m ]- A straightforward calculation yields 



W+[v,m] 



S7T 



d xti 



#(m 1; m 2 )ml 



'm\ + m 2 2 )E{m x , m 2 ) - 2 



m 2 ) 2 



m 



2\2 



where 



H(x,y) 



\og(x 2 /y 2 ) 
x 2 — y 2 



(62) 



(63) 



The symmetry under the exchange of the ordering labels 1,2 is a direct consequence of the 
cyclic property of the trace. 

As noted, this kind of formulas can be evaluated using basis of eigenvectors of m 2 at each 
point x. This gives (Convention 1 is still at work) 



(m 2 + m 2 )H(m v m k 



H(mj, m k )(D tl m RL ) jk (D tl m IjR ) kj 
2 



m 2 k ) 2 



i D » m R)jk{.D ll m Ri 



(64) 



In this formula (D^mRi^jk = {j, L) , etc. The quantities m\ L were defined in 

eq. (D. 

Particular attention requires the case of m\ and m 2 with the same eigenvalue. Before 
momentum integration it is clear that this case is a perfectly regular one. It follows that 
the correct result after momentum integration can be obtained by taking the finite formal 
limit m 2 . — > mf, and then replacing m 2 by its eigenvalue. That is, in eq. (Bl) the case j = k 



(or more generally, m - = m 2 ,) is resolved by taking the limit m 
and finite, namely 



ni 



2 which is well defined 



— / d 2 xT (- 
8vr J 



2 (D^m RL )l 



(65) 



This is the full result in the Abelian case. 

As an application, consider m LR (x) = M(x)U(x), m R L(x) = M{x)U~ 1 {x), where M(x) 
is a c-number but not necessarily constant. Since m 2 = M 2 is c- number, the calculation can 
be done directly from eq. (|62| ) by taking m 2 = M 2 everywhere. This gives 



W. 



2:2 



1 

47T 









j d 2 xtr 











(66) 
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B. Commutator expansions 



The derivative expansion treats m non perturbatively. In other approaches the result 
appears instead in the form of a 1/ m-like expansion such that each term is an homogeneous 
function of m |TT|JT[1. The order of each term in such an expansion is counted by the number 
of commutators (including those implied by the covariant derivative). Therefore, when 
operators m in the middle of an expression are moved to left, producing commutators, the 
new terms so generated count as higher order than the original one and hence they are 
suppressed in this counting. Although the present work is devoted to derivative expansions, 
in this subsection we will consider commutator expansion of our formulas. 

Using our Convention 3 it is remarkably simple to systematically reexpand W + in terms 
of commutators. This can be illustrated with W 22 in form given in eq. flSOl) (this expression 
will be derived below). From eq. ([2"ID it follows that the quantity 



Ci = m x - m 2 (67) 

is equivalent to a commutator [m, ] on the first m^. For instance, 

cjml = (mi - m 2 ) 2 ml = (mi - m 2 )[m, m^m^ = [m, [m, m^m^ . (68) 

Therefore, a commutator expansion is obtained by removing the variable m 2 in favor of c x , 
i.e., m 2 — mi — Ci, and carrying out an expansion in powers of q: 

m 1 m 2 H(m 1 , m 2 ) - 1 = 1 1 1 c t 2 cj 
(mi — m 2 ) 2 6 m\ 6 m\ 15 m\ 

This immediately translates into 

W+ 2 [v,m] =—J d 2 xtr [-L/nJ + -L[m, m^m, + j^[m, [m, mJK + • • 



(70) 



For W± (or higher orders) the identity c n = m n — m n+ i can be used, where c n denotes 
the commutator [m, ] on the n-th fixed element. 

Due to the cyclic property (and in particular integration by parts when is involved) 
neither the derivative expansion nor the commutator expansion take a unique form. How- 
ever, there is an important difference between both expansions regarding their uniqueness, 
namely, the derivative expansion is related to (gauge covariant) dilatations of the external 
fields v and m; each order is tied to a given power of the dilatation parameter. This guar- 
antees that each order in the derivative expansion is a well-defined functional of v and m. 
(That is, each given order can be written in different forms but has a single numerical value 
for each given configuration of the fields.) On the other hand, the commutator expansion is 
not tied to any expansion parameter and so it is possible to have two different commutators 
expansions of a single functional which differ numerically at every order. For instance, in 
the identity 

tr(2m[m, m^m^) = tr([m, [m, m^m^) (71) 

the same functional is represented by two terms of different order (i.e., with a different 
number of commutators), even imposing a standard form (mi and Ci have been chosen 
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as the independent variables). More generally, a functional of the form tr(/(mi, m^mV) 
depends only on the symmetric component of the function f{x\,x 2 ), but the antisymmetric 
component may have a non vanishing commutator expansion when mi and Ci are used as 
the independent variables. This suggest to first symmetrize / and then expand.f] In any 
case, it is clear the usefulness for our notation to analyze this kind of problems and even 
more so for more complicated cases, such as contributions to W± or higher orders. 

The ambiguity in the commutator expansion allows to reorder W 22 so that the symmetry 
under the exchange of ordering labels 1,2 in the right-hand side of eq. (|69|) is restored. A 
convenient choice, containing even orders only, is 

m 1 m 2 -ff(m 1 , m 2 ) — 1 11 \ C\ c 2 1 c\ 



[mi — m 2 ) 2 6 m\m 2 30 m\ m\ 140 m\ m\ 



+ ■••, (72) 



where c 2 = m 2 — m% = — c± represents [m, ] on the second m^. This expansion can be 
obtained from that in eq. ( |6~9|) by recursively subtracting the symmetric version of the leading 
term and then symmetrizing the remainder. Eq. ( |72|) translates into 

W£ 2 [v,m] = — I d-.rU- 



1/1 \ 2 1/1 \ 1/1 \ 

- — m u H — -\m, mj H — -\m, \m, mj + 

6 Vm V 30Vm 2L ' M V 140Vm 3L ' L ' ^ J 7 



(73) 



The commutator expansion just derived is well suited for the vector-like case. In many 
cases, however, an expansion involving commutators of m 2 will be more useful. For instance, 
in the Abelian case m 2 is always a (not necessarily constant) c-number whereas m is not (i.e. 
[m, X] needs not vanish even when X is a multiplicative operator). The method to obtain 



such an expansion is the same as above and can be applied to W 2 2 in eq. (|62|). Namely, the 
quantity 

d = ml - m\ (74) 

represents [m 2 , ] applied to the first fixed factor (m M or (m 2 )^) thus m\ can eliminated in 
favor of C\ in eq. (|62"D. This gives 

^ (mi ' m2) = if + 2^ + 3^ + 4^ + 5^ + ---' 

-2 0/ -3 r-/-4 



{m 2 + m 2 )H(m 1 ,m 2 )-2_ 1 d 3q 2C\ 5Q 



m 2 - m 2 ) 2 ~ 6ml + 6mf + 20m? + 15mi° + 42m} 2 + " ' ' (?5) 



Further, applying a (in this case partial) symmetrization in the labels 1,2 to remove terms 
with an odd number of commutators, yields 

^ 1 fi 1 1 C l C l \ 

H(m u m 2 = — 1 + -^-^ + ^^-^ + ■■■] , 
m{ \ 6 m\m 2 30 m\m\ ) 

(m 2 + m 2 )H(m 1 ,m 2 )-2 1 1 1 C X C 2 1 C\C\ 

+ an ™4„4 ' /ion ~6 „6 + ' ' ' • I '"J 



[mj - ml) 2 6 m\m 2 2 60 m\m\ 420 m\m% 



7 What seems to be true is that, for a given non vanishing functional, the order of the leading 
term (i.e., the number of commutators in the term with the least number of commutators) has an 
upper bound over the set of all possible commutators expansions of the functional. 
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(Where C 2 = m 2 — m| = — C\.) This translates into 



W£ 2 [v,m] 



in 



d xtr 



m: 



6 Vm 2 



30 Vm 4 



■[m 2 ,[m 2 ,m^]]) + 



1 / 1 



6 \m^ m ^ 



— —\m z ,(m z )J 
60Vm 4L ' v m 



2 1 / 1 \ 2 
420 Vm eL ' [ ' v ;mJ V 



In fully expanded notation the leading term is (compare with eq. (|65|)) 

1 



^.leading K m ] 



167T 

1 / 1 



d 2 xti 



—Dnm RL b^m LR + ^D^m LR D^m RL 



m 



R 



6 K , 



1 M A 2 

- — 7rD,,m T 



2 ■ 



(78) 



In the example considered above, m 2 = M 2 a (non necessarily constant) c-number, this is 
the full contribution. 



C. Alternative expressions 



The expression in eq. fl62|) can be somewhat simplified by using the identity (valid inside 
the trace) 



m 



2\2 



{m, m M } 2 = (mi + m 2 )(m 2 + m 3 )m 



>x + m 2 ) 2 m 2 l 



(recall that m 3 = mi due to the cyclic property). This yields the formula 

m 1 m 2 H(m 1 , m 2 ) - 1 



1 

47T 



d 2 xtr 



-m, 



(79) 



(80) 



In most instances, in order to use this kind of formulas it should be clear which are the 
chiral labels L, R in each of the factors. The simplest way to determine this is by forcing 
that the ordering labels 1,2 appear only in m 2 , since this object is either RR or LL and 
so does not flip the chirality label. This can always be achieved by splitting the expression 
into components which are even or odd under mi — ► —mi and m 2 — > — m 2 , that is, in an 
expression of the form tr[/(mi, m 2 )m 2 ] 



/(mi, m 2 ) = A(ml, m 2 ) + m 1 m 2 B(m 2 1 , m 2 ) . 



(81) 



(Note that any expression must be even under m — > — m since there should be the same 
number of R and L labels. Therefore no terms of the form miD(m 2 , m 2 ) + m 2 F(m 2 , m 2 ,) can 
appear.) Applying this procedure to eq. (jS0|) yields 



+ 



1 



d xtr 



[m 2 + m 2 )i/(mi, m 2 ) - 2 



47T 

2mlmlH(m 1 , m 2 ) 



. mt 



m 



m 



2\2 



mm 



mf 



mi) 2 



(82) 
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Now, from our conventions it unambiguously follows that in the first term the chiral labels 
of m 2 and m\ as well as of both factors (mm^) are RR. In the second term, m\ is RR, m\ is 
LL, the first factor is RL and the second one is LR. (Recall that under the Convention 
1, each term is identified with its pseudo-parity conjugate.) This allows to take matrix 
elements as in eq. (|64] ) in a straightforward manner. (The original expression eq. (|62| ) was 
already in this form.) The procedure expressed by eq. (0) shows that any expression, such 
as that in eq. ([$!]), can be brought to a standard form, eq. Q82]), which is free of ambiguities, 
i.e., the expression can be expanded undoing the Conventions 1 and 2 in an unambiguous 
way. This procedure immediately extends to more than two variables, as required in W± or 
higher order terms. 

An essential point in this discussion is that, in an expression tr(/(mi, m 2 )m 2 l ), the func- 
tion f(xi,x 2 ) must be regular in the coincidence limits. To further analyze this point, 
consider the contributions to W 22 in eq. (|80|) coming solely from t/ M , i.e. setting <9 M to zero 
(recall that D^ = + t/ M ). Using the identity 



\2 2 



(83) 



which holds inside the trace, this gives 



WX[v,m] 



1 

Aix 



d xtr 



{m x m 2 H{mi, m 2 ) 



Dv 



+ 0{d) . 



(84) 



Because the action depends exclusively on the combination D — d + v and no algebraic 
assumption has been made on v (it can be an arbitrary matrix), it is clear that all the 
effective action can be reconstructed from the case = by means of the replacement 
v — > D everywhere (this a reciprocal of the usual gauging procedure). So in a formal sense 



W. 



2,2L"> 



— / d 2 xtr 

47T 



'm l m 2 H(m ll m 2 ) - 1) D 2 



(85) 



This is meaningful provided that D M (or i/ M prior to the "gauging" ) appears only in commu- 
tators. Under our Convention 3, this can be easily enforced through the formal identity 



D: 



mi - m 2 



(86) 



The minus sign comes because m M = D^m — mD^ = (m 2 — m^D^ for the first D M in eq. (p5|), 
whereas m /t = (m 3 — m 2 )D M = (mi — m 2 )D^ for the second one. This immediately recovers 
the correct expression eq. (|S0|). 

This does not mean however that the formal procedure is always justified. In fact, we 
could have considered, instead of eq. (0), a different formal expression, namely 



D, 



m 



2\2 



mi — m 



2\2 ■ 



(87) 



This can be used in eq. (^) and the result can be rewritten (in the spirit of eq. 

l r ,o . l 



as 



WvJti, m] = / d z xti . ,, 

2 ' 2L 1 Air J (ml 



ml) 



2\2 



H(m 1} m 2 )(m(m 2 ),) 2 -(m% 



(88) 
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This expression is incorrect if taken naively. Because it involves derivatives of m 2 only, it 
would predict, for instance, a vanishing value for W 2 2 when m is on the chiral circle, i.e. when 
m 2 is a constant c-number, or that W 2 + 2 does not depend on the axial field in the Abelian 
case, both predictions being wrong. The reason is that Eq. (|86D introduces a singularity at 
mi = m 2 which is canceled by the numerator mim2H(mi, m 2 ) — 1, but eq. ( |8"TD introduces 
a new singularity at m\ = — m 2 which is not canceled; each of the terms in eq. (|S8|) is not 
separately finite when m\ = m\ and in general they cannot cancel to each other due to their 
different chiral labels. This renders the expression meaningless unless (m 2 )^ is expanded to 
cancel the spurious singularity at mi = —nri2. 

The fact that the formal expression eq. fl85|) can be promoted to a regular expression 
with D M in commutators (cf. eq. (|80|) ) is a manifestation of gauge invariance. It is perhaps 
interesting to note that the gauge invariance of eq. (|85| ) can be checked even without actually 
bringing D M into commutators. Namely, D M appears only in commutators in an expression, 
if and only if, the expression remains invariant under the shift D M — > D M + a M , where a M is 
an arbitrary constant c-number. That the formula eq. fl85|) is gauge invariant can be seen 
by applying this test: it is readily seen that in the terms with a M , /r?i and m 2 coincide (since 
m and a M commute) and the limit mi — > nri2 gives zero. (That it is a double zero follows 
from symmetry under exchange of labels 1,2.) 

Finally, let us mention that eq. (0), and thus eq. (|85p, can be obtained rather directly 
in several ways without using Chan's method. The key point is that in the absence of <9 M 
all quantities are multiplicative operators and the cyclic property holds. The procedure 
implied by combining eq. ( p5|) and eq. (|86|) , can then be extended to compute H^t^m] i n 
any number of dimensions. However this method is not suitable to treat the case of four or 
more derivatives. For instance, the formula similar to eq. ( |35l) to fourth order is 

w 4,2 = 4- J d 2 xtr ((mim 3 - m 2 m 4 )iy 4i i + m 1 m2m 3 m i H 4fi ) D^ t D fl D u D u 

~ ((mi/113 + m 2 m A )H 4A + mim 2 m3m4# 4i o) D^D U D^D V , (89) 

where are the functions of m 1; m 2 , m 3 , m 4 defined in eq. (^0]). Using their explicit form 
and the test noted above of shifting D M by a constant c-number, it is possible to check that 
the previous formula is chiral invariant (of course, this is just a check of the calculation). 
However, there is no simple systematic form for bringing it to an explicitly invariant form 
where all covariant derivatives appear only in commutators. The trick introduced in eq. flSSP 
works in W 2 + because appears there as a first covariant derivative. Such replacement is 
no longer sufficient in since explicit chiral invariance requires at least the presence of 
second covariant derivatives as well as F (If only the replacement in eq. (|86|) is used, one 
obtains an expression involving with coefficient which are functions of mi, m 2 , m±. 
However, these functions are not finite in the coincidence limit, i.e. when two or more of 
their argument become equal, and such an expression is not truly well-defined.) 



D. Bosonic formulas 



For completeness we give here the analogous formulas of eqs. (|3^j38|j3^) and eqs. (|59|j60|j6T|) 
for the bosonic effective action. 
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m , / ^9 f d d xd d p 



(27r) d 
2p 4 



d(d + 2) [ 

/ 1 ^ 2 
+2 ^ 



p 2 / 1 ^ 2 

/ 1 \ 4 / 1 1 ^ 2 



1 1 ^ 2 



, 1 1 1 



1 1 1 / 1 \ 2 1 



(90) 



where A = p 2 + U, U il = [P M , [/], etc. From this formula, the fermionic one is obtained by 
setting U = m 2 + + . 

After momentum integration, the renormalized bosonic effective action becomes 



Trlog(P 2 + f/) 



f 



(4n) d / 2 T(d/2) 
4 



/ d d X t,V-lf 2 \u d l 2 \o g (U/^ 2 ) + \H% /2 Ul 



H} 2 J' 2 M 2 U 2 



— H 2 ' 2 ' 2 ' 2 (U II V 

2) -"4,(2/2+1 l^M^ J 



d(d + 2) 

I() -3,1,3,1 rr 2 rr 2 4 u-3,3 rr 2 16 H^lt , -JJ 2 U V 



rr3, 1,3,1 tt2tt2 JJ^,3 TT^ _i_ rr3,l,3 

"d(d + 2) ^V/2+l^ rf(d + 2)^2,^/2+1^ + d{d+ 2) ^3,d/2+ 



Tt2,2,2 7-7- 7-7- p 

- n 3,d/2+l u M u i'- r Mf 



-if. 



2,2 



-l r (* 



d(d + 2) " W+i - "~ ^ d(d + 2) JJ 2 ><W ^ ^ 
In this case, the arguments z 2 , . . . , z 2 in H^ 1, k "' Vn have to be substituted by C/i, . . . , U n . 



(91) 
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